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Context. The classification scheme for high redshift galaxies is complex at the present time, with simple colour selection criteria 
f*"^ ' (i.e. EROs, IEROs, LBGs, DRGs, BzKs), resulting in ill-defined properties regarding stellar mass and star formation rate for 
t . these distant galaxies. 

'wJ ' Aims. The goal of this work is to investigate the properties of different classes of high-z galaxies, focusing in particular on the 
Stellar Masses of LBGs, DRGs and BzKs, in order to derive their contribution to the total mass budget of the distant Universe. 
Q ■ Methods. We have used the GOODS-MUSIC catalog, containing ~3000 Ks-selected (~10000 z-selected) galaxies with multi- 
wavelength coverage extending from the U band to the Spitzer 8[im band, with spectroscopic or accurate photometric redshifts. 
We have selected samples of BM/BX/LBGs, DRGs and BzK galaxies, discussed the overlap and the limitations of these criteria, 
which can be overcome with a selection criterion based on physical parameters. We have then measured the stellar masses of 
these galaxies and computed the Stellar Mass Density (SMD) for the different samples up to redshift ~ 4. 

Results. We show that the BzK-PE criterion is not optimal to select early type galaxies at the faint end. BzK-SF, on the other 
hand, is highly contaminated by passively evolving galaxies at red z — Ks colours. We find that LBGs and DRGs contribute 
k>( \ almost equally to the global SMD at z > 2 and in general that star forming galaxies form a substantial fraction of the universal 
j_j . SMD. Passively evolving galaxies show a strong negative density evolution from redshift 2 to 3, indicating that we are witnessing 
03 the epoch of mass assembly of such objects. Finally we find indications that pushing the selection to deeper magnitudes, the 
contribution of less massive DRGs could overcome that of LBGs. Deeper surveys, like the HUDF, are required to confirm this 
suggestion. 
Conclusions. 

Key words. Galaxies:distances and redshift - Galaxies: evolution - Galaxies: high redshift - Galaxies: fundamental parameters - 
Galaxies: mass function 
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1. Introduction 

Until recently, the Lyman drop-out technique was the only 
selection criterion used to identify large numbers of galax- 
ies at z > 3 in deep optical surveys (jMadau et al. 1 996, 
Steidcl ct al. 2003}. This technique is effective in find- 
ing star forming galaxies, characterised by a typical blue 
spectrum with reduced dust absorption, the so called 
Lyman Break Galaxies (LBGs). In recent years, with the 
availability of large IR detectors at 4-8 meter class tele- 
scopes, the ISO and Spitzer data in the medium IR and 
the SCUBA observations, new types of high-z galaxies 
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have been identified, with spectral energy distributions 
(SED) and physical properties (Age, Mass, SFR, dust 
content) different from the "classical" LBG population. 
Among them, it is worth reminding the sub-mm galax- 
ies (Sm ail et al. 2004|) . the distant red galaxies (DRGs, 
IFranx et al. 20031 Ivan Dokkum et al. 2003]) . Extremely 
Red Objects (EROs, |McCarthy 20041 IDaddi et al."2000l) , 
IRAC Extremely Red O bjects (IEROs, lYan et al. 2004|) 
and the "BzK" galaxies (|Daddi et al. 2004|) . The aim of 
these different selection criteria is to find objects that 
could be missed by the Lyman Break technique, due to 
the presence of an evolved stellar population (dead and 
massive objects) or a young stellar population heavily ob- 
scured by dust at high redshift (dusty starbursts). 
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While the LBG selection identifies galaxies according 
to their less obscured star formation properties based on 
UV rest frame colours, EROs, DRGs, IEROs and BzK 
selection criteria include also optical rest frame colours, 
which are less affected by recent SF episodes, but are much 
more sensitive to the light coming from lower-mass longer- 
lived stars. Another advantage of the selection from the 
near-IR is that the rest frame optical light is less affected 
by dust extinction. However it is not yet clear how these 
different types of galaxies overlap and what are the galax- 
ies that dominate the number counts and the mass density 
budget at high redshift. 

Using different SFR indicators, |Reddy et al. 2005| and 
|Reddy et al. 2006] discussed the properties of LBGs, 
DRGs and BzK galaxies in the redshift range 1.4 < z < 
2.6 from the GOODS-North survey, and found a strong 
overlap (70-80%) between LBG and BzK galaxies, when 
the analysis is restricted to a Ks-selected sample. 

LBGs, DRGs and BzKs present a similar SFR distri- 
bution as a function of the Ks magnitude, the M/L ratio 
of DRG and BzK galaxies (both massive and dusty) are 
larger than LBGs, but their stellar mass does not exceed 
the range spanned by optically selected galaxies. Previous 
studies ( |Adelberger et al. 200"4"l|Shapley et al. 2005] ) point 
to a similarity in the metallicities, clustering and stellar 
masses of Ks < 20(vega) optical and near-IR selected 
galaxies at 1.5 < z < 2.5. Indeed, Red dy et al. 2 005 con- 
cluded that the presence or absence of star formation may 
be the only significant difference between optical and near- 
IR selected massive galaxies (M star > 1O 11 M ), and the 
difference in SFR may be temporal, i.e. due to a phe- 
nomenon of transient SF activity. On the other hand, 
IQuadri et al. 2006] have studied the clustering properties 
of K-selected galaxies at 2 < z < 3.5 in the MUSYC sur- 
vey and found that their correlation length does not de- 
pend on the K band magnitude, but it increases strongly 
both with J — K and with R — K colours. This suggests 
that K-bright blue galaxies and K-bright red galaxies are 
fundamentally different, at least in the clustering proper- 
ties and hence on their host halo mass. 

Recently, Ivan Dokkum et al. 20061 compared the DRG 
and LBG selection criteria at 2 < z < 3 using a K- 
selected sample (Ks < 21.3(vega) corresponding to M > 
1O U M ), finding that DRGs make up 69% (77%) of the 
sample by number (by mass) and LBGs only 20% (17%), 
with a small overlap between the two groups. The cen- 
sus of the stellar mass density at high redshifts and faint 
luminosities however is not well sampled in the recent lit- 
erature as the SFR is, probably due to the lack of a com- 
prehensive database for LBGs, DRGs and BzKs, where ac- 
curate redshift measurements and stellar mass estimates 
are available. 

The GOODS-MUSIC sample (|Grazian et al. 2006aj) 
allows to investigate the role of different galaxy types 
(LBGs, DRGs and BzKs) in the universal stellar mass 
budget at high redshift. Using this sample we are 
able to select on the same wide, deep area all galaxy 
types, verify the effectiveness of different colour crite- 



ria in sampling galaxies at high redshifts, study the 
overlap within the sub-samples and their contribution 
to the total stellar mass density of the Universe. In 
this field, the stellar mass function for all galaxy types 
are already studied QDrory et al. 2005 Caputi et al. 2 004 



Font ana et al. 2006]) . the latter taking advantage of 
the IRAC-Spitzer data, which are fundamental to de- 
rive the stellar mass at higher-z ( |Papovich et al. 2 005, 
Fontana et al. 2006]) . We focus here only on two simple 
but largely debated questions, i.e. the role of LBGs, DRGs 
and BzKs to the total stellar mass density of the Universe 
at z > 1.4, and the comparison between optically selected 
and near-IR selected samples. 

The paper is organised as follows. In Sect. 2, we remind 
the basic feature of our GOODS-MUSIC dataset. In sec- 
tion 3 we describe the selection criteria used to identify 
LBGs, DRGs and BzK galaxies, the overlap between the 
samples, the observational limitations and their intrinsic 
properties. In Section 4, we present the basic results of 
our analysis, namely the resulting mass density per galaxy 
type, and its redshift evolution. In Section 5 we discuss the 
results and derive our conclusions. 

All magnitudes are in the AB system (except where 
otherwise stated) and we adopt the A-CDM concordance 
cosmological model (Hq = 70, Qm = 0.3 and f^A = 0.7). 



2. The Data 

In this work we use the GOODS-MUSIC (GOODS 
Multicolour Southern Infrared Catalog) sample, a 14 
bands multicolour catalog extracted from the deep and 
wide survey conducted over the Chandra Deep Field South 
(CDFS, Giavalis co et al. 2 004). in the framework of the 
GOODS public survey. The data comprise a combination 
of images that extend from U to 8.0/im, namely U-band 
data from the 2.2ESO {U 35 and U 3S ) and VLT-VIMOS 
(Uvimos), the F435IF, F606W, F775W and F850LP 
(Z-band) ACS images, the JHKs VLT data and the 
Spitzer data provided by the IRAC instrument at 3.6, 4.5, 
5.8 and 8.0 //m. The total area of the GOODS-MUSIC 
database results from the overlap of ACS F435W band 
and Ks ISAAC images, for a total of 143.2 sq. arcmin., 
which is covered by 12 bands, most notably the Spitzer 
ones (U 35 U 38 BVizJKs and IRAC). Only the U-VIMOS 
and H bands have limited extension, of 90.2 and 78.0 sq. 
arcmin., respectively. 

Since the detection mosaics have a complex, inhomoge- 
neous depth and are quite different in the z and in the Ks 
bands, we have divided the whole sample into 6 indepen- 
dent catalogs for the z band and in 6 for the Ks band, each 
with a well defined magnitude limit and area, that we use 
to compute the stellar mass densities and other statistical 
properties in this paper. The limiting magnitudes for each 
area are reported in Table 2 of Grazian et al. 2006a and 
refer to a 90% completeness. The typical magnitude limit 
for most of the sample is about Ks = 23.5 and z = 26.0 
(corresponding to an area of 72 and 99 sq. arcmin., respec- 
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tively), and extends down to Ks — 23.8 and z = 26.18 in 
limited areas, as described in lGrazian et al. 2006al 

The GOODS-MUSIC database contains both a z- 
selected and a K s-selected catalog. Colours have been 
measured using a specific software for the accurate "PSF- 
matching" of space and ground based images of different 
resolution and depth, that we have named ConvPhot 
(De Santis et al. 2006). We have cross correlated our 
catalog with the whole spectroscopic catalogs available 
to date, from a list of surveys, assigning a spectroscopic 
redshift to more than 1000 sources. In this work we use 
a spectroscopic sample that is wider than that presented 
in IGrazian et al. 2006a|, thanks to the increased num- 
ber of spectra publicly available (Vanzcll a et al. 2006|) . 
Finally, we have applied our photometric redshift 
code, developed and tested over the years in a series 
of works (jFontana et al. 20001 ICimatti et al. 2002| 

IFontana et al. 2003| IFontana et al. 2004| 

|Giallongo et al. 2005[ IFontana et al. 2006|) . that adopts a 
standard x 2 minimisation over a large set of templates 
obtained from synthetic spectral models. The comparison 
with the spectroscopic sample (jGrazian et al. 20 06a, 
IGrazian et al. 2006bp shows that the quality of the 
resulting photometric redshifts is excellent, with a r.m.s. 
scatter in Az/(1 + z) of 0.03 and no systematic offset 
over the whole redshift range < z < 6. The procedures 
that we adopted to extract this catalog and derive 
the photometric redshifts are described at length in 
IGrazian et al. 2006a 

The final samples that we adopt here consist of 2931 
galaxies down to Ks ~ 23.8 and 9862 galaxies down to 
z = 26.18, taking into account the complex magnitude 
limits of the survey. On the total if s-selected sample, 
815 galaxies (973 for the z-selected) have reliable spec- 
troscopic redshifts and the remaining fraction have well 
trained photometric redshifts. 

AGNs and stars are identified first from spectroscopic 
information. We then distinguish galaxies from stars and 
AGNs using morphological and photometric information, 
when spectroscopic data are not available. Point-like 
sources are selected using the star/galaxy separation flag 
(s/g) provided by SExtractor (jBertin &: Arnouts 1996ft in 
the z band. We tune the selection on known spectroscopic 
stars, as described in lGrazian et al. 2006al We use photo- 
metric information to further check this criterion, using in 
particular the "BzK" colour criteria of Da ddi et al. 2 004 

For the purposes of the present work, we will use the 
-fTs-selected catalog to study the properties of DRGs and 
BzK galaxies, while the LBG sample is derived from the z- 
selected catalog. When discussing the properties of LBGs, 
we restrict the GOODS-MUSIC sample only to the area 
covered by deep U-VIMOS observations, since LBGs are 
selected by means of U, V and I bands. In practise, for 
LBGs, we use the area of 90.2 sq. arcmin. with U-VIMOS 
coverage, while for DRGs and BzK galaxies, we use the 
total GOODS-MUSIC area (143.2 sq. arcmin.). 
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3. Similarities and differences between LBGs, 
DRGs and BzK galaxies 

We will now compare the various colour criteria used to 
find high-z galaxies. In particular we will discuss the sim- 
ilarities and differences between LBGs, DRGs and BzK 
galaxies, to shed light on two important questions, namely 
the comparison between z-selected and if s-selected sam- 
ples of galaxies at high redshift, and the relative fraction of 
different galaxy types in a K s-selected sample of galaxies. 

We remark here that these simple colour criteria 
(LBGs, DRGs and BzKs) cannot distinguish between pas- 
sively evolving and evolved galaxies. A passively evolving 
galaxy is defined as an object where SF has ceased by 
a time much longer than the typical duration of the SF 
episodes, while an evolved galaxy has a large population 
of old stars but can also have recent /on-going star forma- 
tion. Usually, through simple colour criteria, it is possible 
only to check for the presence of evolved stellar popula- 
tions. If the wavelength baseline is extended, as in the case 
of the GOODS database, it is possible to isolate also true 
passively evolving galaxies. 

3.1. Selection of LBGs, DRGs and BzK galaxies 

An efficient method to select unobscured and modestly 
obscured star forming galaxies at high redshifts is the 
Lyman break technique, that is effective at 2.8 < z < 
3.7, as defined originally by ISteidel fc Hamilton 1 993, 
IMadau et al. 19961 ISteidel et al. 19951 ISteidel et al. 19991 
A recent extensions at lower redshift are the so called 
"BX" galaxies, at 2.2 < z < 2.8 and "BM" galaxies, 
at 1.4 < z < 2.2 flAdelberger et al. 2004[ ). Since the fil- 
ter set of the GOODS-MUSIC sample is different from 
the classical UGR adopted by Steidel and collaborators 
(|Steidel et al. 2 003). we have used the following criteria, 
tuned by a comparison to the photometric redshift distri- 
butions: 

-0.2 < V-I < 0.35 
U - V > 0.75(V - 1) + 1.15 (1) 

for LBGs, 

£7 - F > 0.65(F - 1) + 0.25 
U — V < 0.75(V — I) + 1-15 

-0.2 < V-I < 0.25 (2) 

for "BX" and 

u -v >om(v -i) 

[/ - V < 0.65(U - /) + 0.25 

-0.2 < V-I < 0.45 (3) 

for "BM". 
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We have also verified that these criteria are almost 
equivalent to the original BM/BX/LBG criteria by repro- 
ducing for each galaxy in our catalog the synthetic U GR 
magnitudes of Steidel and collaborators using the best fit 
spectrum derived by our SED fitting technique. The colour 
selections used here are also consistent with criteria based 
on photometric redshifts. 

The differences in the filter sets, colour selec- 
tions and the fact that the BM/BX/LBGs are tra- 
ditionally selected in the R band and not in the 
z band, as in this work, result in a rcdshift distri- 
bution of the BM/BX/LBG galaxies wider than in 
Adelberger et al. 2004| and |Reddy et al. 2005[ There is a 
small sample of BX and LBGs at z p hot ~ 0.4, due to galax- 
ies at lower redshifts with U — V and V — I colours typi- 
cal of high-z galaxies. In the following, however, we study 
the properties of BM/BX/LBGs selected with our colour 



criteria only at 



Zphot 



> 1.4, thus avoiding the contam- 



ination from these lower redshifts galaxies. In addition, 
BM/BX/LBGs were originally relatively bright galaxies, 
since the historical criterion used was R < 25.5, while in 
this work the selection is pushed till the nominal com- 
pleteness of the GOODS-MUSIC sample. Adopting for 
our BM/BX/LBG sample a fainter limit in a redder band, 
however, should not introduce further contaminants, since 
these galaxies are selected for their typical blue optical 
colours. 

When applied to the i^s-selected sample, these colour 
cut, with the additional constrain z > 1.4 for consistency 
with the classical BM/BX/LBG criteria, selected 166 
galaxies with substantial star formation at 1.4 < z < 3.7 
and Ks < 23.8, with U-VIMOS magnitude determination: 
hereafter we will refer to this sample as BM/BX/LBG-Ks. 
The z band, however, is much more appropriate to select 
star forming galaxies (with low dust extinction), though 
it is not the best choice for the stellar mass estimate, in 
particular at high redshifts. Using the same criteria as 
above, we have therefore selected 1345 galaxies in the z 
band (z < 26.18), and we will refer to this sample as 
BM/BX/LBG-Z. 

Recently, IDaddi et al. 20041 proposed a criterion to 
identify galaxies at 1.4 < z < 2.5, differentiating them 
between passively evolving (BzK-PE), selected according 
to the criterion BzK = (z — K)ab — (B — z)ab < —0.2 
and (z-K)ab > 2.5, and star-forming galaxies (BzK-SF), 
which can be effectively isolated from lower redshift inter- 
lopers and stars through the simple criterion BzK > —0.2. 
The advantage of the BzK criterion is that this colour 
combination is insensitive to extinction since the redden- 
ing vector is parallel to the BzK = —0.2 dividing line, as 
shown in lDaddi et al. 20041 This technique reveals partic- 
ularly efficient for sampling the dusty side of star forming 
galaxies at intermediate/high redshifts. Using the BzK se- 
lection criteria, we have extracted 89 BzK-PE and 747 
BzK-SF galaxies down to the conservative Ks complete- 
ness limits of the survey, described in Section 2 and in 
IGrazian et al. 2 006a 
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Fig. 1. An example of LBG, DRG, BzK-SF and BzK-PE 
from the GOODS-MUSIC sample (jGrazian et al. 2006ap . 
The x-axis corresponds to the observed wavelengths, while 
the y-axis refers to the observed magnitudes in the AB 
photometric system. 



Finally we focus on the so-called Distant Red Galaxies 
(DRGs, IFranx et al. 2 003) . These galaxies are selected 
through the J — K > 2.3(vega) criteria, originally 
designed to be sensitive to galaxies with large 4000 
A and/or Balmer breaks at z > 2. These features 
become strong for ages larger than 1 Gyr, even if 
they are also sensitive to the metallicity of galaxies 
dPoggianti fc Barbaro 19971 IKaufi mann et al. 2002] ). The 
presence of a Balmer and/or 4000 A break can thus 
be an indication of an evolved stellar population. It 
turned out also that this single colour selection is sen- 
sitive both to evolved galaxies and to dusty starbursts 
at z > 1 — 2, as shown in Forster Schreib er et al. 2 004, 
|Papovich et al. 2005| and IGrazian et al. 2006bl DRGs 
thus appear to be a mixed/composite sample. The DRGs 
in the GOODS-South sample are already well studied 
(Papovich et al. 2005 Grazian et al. 2006b): we select 179 
galaxies of this type down to Ks = 23.8, which span a 
large redshift range 1 < z < 4. 

In Fig. Q] we show typical SED of an LBG, a DRG, a 
BzK-SF and a BzK-PE galaxy. 



3.1.1. The number density of LBGs, DRGs and BzK 
galaxies in the GOODS South field 

The area of the GOODS South field is relatively large 
compared to the deep pencil beam surveys like the two 
Hubble Deep Fields and the Hubble Ultra Deep Field, 
but it is still limited and therefore subject to cosmic vari- 
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ance, as demonstrated by the prominent peaks in the spec- 
troscopic redshift distribution (Vanzcll a et al. 2006)) . It is 
also known that the number density of X-ray sources in 
the GOODS-South is lower compared to other deep X-ray 
surveys, i.e. the Chandra Deep Field North at faint fluxes 



and hard energy band (2-8 keV, |Brandt fc Hasinger 2005 ), 
and there are claims that this underdensity extends to op- 
tically/NIR normal galaxies. 

We have derived the number densities of 
BM/BX/LBGs, DRGs and BzK galaxies in the GOODS 
South and compared them to the densities obtained by 
large (but shallower) area surveys, in order to verify if 
GOODS-MUSIC is a representative sample. 

In their 0.81 sq. deg. survey, Adelberger et al. 2004 
found that at R < 25.5 the number densities of BM, BX 
and LBGs are 3.28, 4.82 and 1.70 galaxies per sq. arcmin., 
respectively. For comparison, in the GOODS-MUSIC sam- 
ple our criteria for BM, BX and LBGs select 3.35, 5.64 and 
2.22 galaxies per sq. arcmin., at a limiting magnitude of 
Z < 25.0, Z < 25.3 and Z < 25.33, respectively. If we limit 
our BM/BX/LBGs at z p hot > 0.5 the number densities are 
3.22, 4.40 and 1.70 galaxies per sq. arcmin. and therefore 
agree much better with those of |Adelberger et al. 2004| 

Recently, [Quadri et al. 2006] found 0.89 DRGs per sq. 
arcmin. at Ks < 2\{vega); in the GOODS-MUSIC sam- 
ple at the same magnitude limit the number density 
is similar, with 0.80 DRGs per sq. arcmin. Indeed, the 
DRG counts in the GOODS South is consistent with 
that derived in deeper surveys, like FIRES, as shown in 
IGrazian et al. 2006bl 

The seminal work for the BzK selection criteria by 
Daddi ct al. 2004 is based on the K20 survey (a subset 
of the GOODS-South field) and the number densities of 
BzK-SF and BzK-PE at K < 20(vega) are 0.91 and 0.22 
galaxies per sq. arcmin., respectively. In a larger field, 
|Kong et al. 2 006 found, at the same magnitude limit, 1.20 



BzK-SF/c 



and 0.38 B zK — PE j arcmin , respec- 



tively In the GOODS-MUSIC sample, at the same magni- 
tude limit, the BzK-SF galaxies have a slightly lower num- 
ber density (0.77), while the BzK-PE counts are consistent 
(0.39) with the larger area survey of |Kong et al. 2 006. 

An independent estimate of the number densities for 
BzK galaxies is provided by Re ddy et al. 2006} in the 
GOODS-North the number densities are 3.1 and 0.24 
arcmin~ 2 for BzK-SF and BzK-PE, respectively, at K < 
21(vega). In the GOODS-MUSIC sample at this deeper 
magnitude cut, the densities are comparable or slightly 
higher, with 3.2 and 0.65 galaxies per sq. arcmin. for BzK- 
SF and BzK-PE, respectively. The large variation in num- 
ber density for the BzK-PE sample is due probably by the 
large angular clustering of this evolved population and the 
limited area which is more sensitive to the field-to-field 
variations. 

We thus confirm that the GOODS-South region is not 
underdense at intermediate/high redshifts (1 < z < 4) 
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1 The R-Z colour of BM/BX/LBGs changes with redshifts, 
galaxies becoming redder at lower redshifts. 



(B-z) 

Fig. 2. The B — z vs z — Ks colour for K s-selected 
galaxies in the GOODS-MUSIC sample (small dots). We 
do not plot AGNs or stars. Blue small points show the 
star forming galaxies selected with the BM/BX method 
(from the z -selected sample), magenta empty squares are 
the Ks-selected BM-BX galaxies, while red stars are the 
DRGs. Horizontal arrows indicate that the B magnitude 
estimation is based on an upper limit at la. The dashed 
line shows the separation between galaxies (up) and stars 
(down) according to lDaddi et al. 20041 The diagonal solid 
line at BzK — — 0.2 allows to isolate star forming galax- 
ies at 1.4 < z < 2.5 (up) from galaxies at lower redshifts 
or from passively evolving one in the same redshift range 
(z - Ks > 2.5). The reddening vector for E(B - V) = 0.3 
is plotted at the left of the diagram, and is parallel to the 
BzK = —0.2 line. This ensures that the BzK-SF criterion 
is insensitive to dust extinction. 



and, though affected by cosmic variance, it is representa- 
tive of the general properties of the distant Universe. 

3.2. The overlap between samples 

Fig. [H shows the distribution of BM/BXs and DRGs in 
the B — z vs z — Ks diagram. In the BzK plot the loci 
of star-forming (BzK-SF) and passively evolving galaxies 
(BzK-PE) are separated by the BzK = —0.2 diagonal 
line. 

In Table Q] we report the overlap between different 
populations: these samples are restricted to the area cov- 
ered by U-VIMOS and limited for if s-selected galaxy at 
1.4 < z < 2.5, in order to compare homogeneous classes. 
We add to TabfT]the BM/BX galaxies that are ^-selected, 
for the sake of completeness. 
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Table 1. Overlap between DRG, BM/BX and BzK galaxies at 1.4 < z < 2.5 



Sample 


DRG 


BzK-PE 


BzK-SF 


BM/BX-Ks 


BM/BX-Z 


DRG 


78 


5/78 (6%) 


72/78 (92%) 


2/78 (3%) 


1/78 (1%) 


BzK-PE 


5/37 (14%) 


37 


0/37 (0%) 


0/37 (0%) 


0/37 (0%) 


BzK-SF 


72/327 (22%) 


0/327 (0%) 


327 


123/327 (38%) 


121/327 (37%) 


BM/BX-Ks 


2/124 (2%) 


0/124 (0%) 


123/124 (99%) 


124 


122/124 (98%) 


BM/BX-Z 


1/913 (0.1%) 


0/913 (0%) 


121/913 (13%) 


122/913 (13%) 


913 



All the samples are restricted to the area covered by U-VIMOS in order to select BM and BX galaxies and limited for all 
galaxy type to 1.4 < z < 2.5 where the BzK criterion is efficient. The BM/BX-Ks sample represents galaxies in the Ks-selected 
database which resemble BM or BX colour criteria, while BM/BX-Z refers to the z-selected sample. The diagonal of the matrix 
indicates the total number of galaxies selected for each type. It is useful to remind that the BM/BX-Z sample is z-selected, and 
thus cannot be compared directly with the Ks-selected samples, but it is useful only as reference. 



From Figure and Table [T] we can draw the following 
considerations. 

Clearly, DRGs and BM/BXs are quite orthogonal as 
selection criteria, with a very small overlap of only 2 (out 
of 78 DRGs and 124 BM/BXs) galaxies in common be- 
tween them. 

On the other hand there is significant overlap between 
BzK and the BM/BX or DRG criteria: 

— The BzK-SF has overlap both to DRGs ( presum- 
ably their subset of dusty star forming galaxies) and 
BM/BXs ( less obscured star forming galaxies): we can 
therefore conclude that it is a good selection criterion 
for star forming galaxies, and quite complete, regard- 
less of dust obscuration properties. 

The fraction of BM/BX galaxies which are also BzK- 
SF is 99%, while the fraction of BzK-SFs that are also 
BM/BXs is 38%. Thus, BM/BX can be seen as a sub- 
sample (less obscured) of typical starforming galaxies, 
represented by the BzK-SF class. This result is slightly 
comparable to that of the GOODS-North field where 
|Reddy et al. 2005| found an overlap of 60 - 80% be- 
tween BM/BX galaxies and BzK-SF at z ~ 2, based 
on a if-selected sample: the overlapping fraction also 
increases towards fainter galaxies. 

— The DRG population has overlap with all other 
selection criteria, with a predominance of BzK- 
SF. Indeed, DRGs are a mix of two populations, 
where star-forming galaxies strongly obscured by 
dust and massive/evolved galaxies coexist, and pos- 
sibly share both these properties. It is well known 
dPapovich et al. 2005| IGrazian et al. 2006bp that low- 
z DRGs are mainly galaxies caught in their dusty star- 
burst phase, which have thus distinct SF properties 
in respect to passively evolving galaxies. The high-z 
DRG subset, instead, is a mix of dusty starburst and 
evolved massive galaxies, as shown by recent deep IR 
spectroscopy (jKriek et al. 2006|) . 

Note however that a number of DRG galaxies that lay 
in the upper left part of the BzK plot, only have limits 
in the B-z colour. Therefore it remains unclear whether 
they would still lay in the BzK-SF locus or migrate to 
the BzK-PE area if deeper B-band data were added. 



In the next section we will discuss these observational 
limitations in more detail. 
— Among the BzK-PE galaxies, 14% are also DRGs, as 
one would expect that evolved galaxies are also se- 
lected by the Balmer and/or 4000 A break criterion. 
However, a lot are not DRGs: from an inspection of 
their SED, we find that these objects are all evolved 
galaxies, with a Balmer and/or 4000 A break at red- 
shift smaller than 2, or with a J — K colour just below 
the defining DRG cut. Indeed, also Rcd dy et al. 2 005, 
based on a spectroscopic selected sample, find that the 
redshift distributions of BzK-PE galaxies and DRGs 
have very little overlap. 

The BzK-PE criterion is orthogonal to methods se- 
lecting star forming galaxies, thus it is not contami- 
nated by obscured starbursts mimicking the colour of 
a red/evolved galaxy. It is therefore an effective crite- 
rion to find evolved galaxies at z ~ 2, although as we 
will discuss in the next section, it is highly incomplete 
mainly due to observational limitations. 

3.3. Observational limitations of LBG, DRG and BzK 
criteria 

In Figure [5] there are a few star forming galaxies selected 
by the BM/BX-Ks criteria which stand below the BzK = 
-0.2 find! This seems to indicate a small incomplete- 
ness of the BzK-SF selection towards star-forming galaxies 
with low extinction, as confirmed by recent VIMOS spec- 
troscopy for a large (> 1000) sample of galaxies in the 
GOODS-South field flPopesso 20061 . However, part of the 
effect is also due to the different redshift interval of the two 
samples, 1.4 < z < 2.5 for the BzK-SF and 1.4 < z < 2.8 
for the BM/BX criterion. 

As noted in the previous Section, most DRGs lay in the 
locus of BzK-SF galaxies, with only a small fraction falling 
in the BzK-PE selection area. However when we analyse 
the galaxies individually we find several DRGs with very 
red z — Ks colour and upper limit in the B band in the 



2 Clearly, the number of z-selected BM/BX galaxies below 
the colour cut is much higher. This is better shown in the colour 
figure where the BM/BX-Z and BM/BX-Ks samples are clearly 
distinguished. 
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locus of BzK-SF, whose SEDs are better fit by an evolved 
galaxy at redshift > 1.4. In agreement with Fig. 12 of 
|Reddy et al. 2005 we find that in the BzK diagram the 



locus proposed bv lDaddi et al. 2 0041 for evolved galaxies is 
not well defined, when considering galaxies with very red 
z — Ks colours. A motivation for this effect is the lack of B 
band deep enough to provide a good diagnostic when very 
red galaxies in z — Ks are considered. Dad di et al. 2 004 
in fact used this photometric criterion on the K20 sample, 
with Ks-band magnitude limited at Ks = 20(vega) and 
taking advantage of the relatively deep B-band image pro- 
vided by GOODS. When fainter galaxies in the Ks-band 
are included, as in this work, the BzK criterion starts to be 
ineffective in isolating early type galaxies at high redshifts, 
as discussed also in |Reddy et al. 20"05 and IRenzini 20061 
In addition, few galaxies with detection in the B band 
and very red z — Ks colour are also better fit by evolved 
galaxy SED, suggesting that there is an intrinsic mix of 
star- forming and evolved galaxies at z — Ks > 2.5. This 
is predicted also bv lDaddi et al. 20041 who in Fig. 8 show 
that galaxies with high formation redshift (zf orm > 5) 
and relatively short e-folding time of the star formation, at 
redshift z > 2 fall in the BzK-SF region before migrating 
to the BzK-PE region at z ~ 1.4. We have checked some 
of our evolved galaxies in the BzK-SF region and they 
have indeed very high formation redshift and z p hot > 2 
(or spectroscopic when available). 

In Fig. [3] we show the photometric redshift distribution 
of BM/BX/LBGs (^-selected), DRGs, BzK-SF and BzK- 
PE galaxies, all limited to the area covered by U-VIMOS 
observations. When available, spectroscopic redshifts are 
used. In this plot there are BzK-SF galaxies at z > 3, 
even if the BzK criterion should be efficient in the range 
1.4 < z < 2.5. A detailed analysis of the SED and physical 
parameters of these objects indicates that they are almost 
all evolved galaxies with the Balmer and/or 4000 A break 
falling between the J and the Ks band, thus resembling 
the DRG criterion. In the BzK diagram they occupy the 
DRG locus with extreme colour in z — Ks. 

The photometric redshift distribution shown in Fig. 
[3] is quite different from that shown by |Reddy et al. 2 005, 
since photometrically selected DRGs and BzK-SFs show a 
wider redshift range than spectroscopically selected ones. 
These differences probably arise because of the limitation 
of the spectroscopic redshifts, which are biased toward op- 
tically bright objects and can depend on the instrumenta- 
tion used to obtain the spectrum and to the spectral type 
of the galaxy. In addition, part of the effect could be due 
to field-to-field variance, which is not negligible even with 
the volumes covered by the GOODS survey. 

In summary, we have carried out a detailed compar- 
ison between the LBG, DRG and BzK selection criteria. 
The selection of galaxies according to the LBG criterion 
is sensitive only to the moderately obscured SF galaxies, 
missing the dusty starburst objects (60% of star forming 
galaxies at 1.4 < z < 2.5). DRGs, instead, are less sen- 
sitive to dust obscuration effects, but comprises a mix of 
two populations, the old/evolved galaxies and the dusty 
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Fig. 3. Photometric redshift distribution of 
BM/BX/LBGs (^-selected), DRGs, BzK-SF and BzK-PE 
galaxies, all limited to the area covered by U-VIMOS 
observations. We limit the BM/BX/LBG-Z sample at 
Zphot > 0.6 to avoid the lower redshift contaminants 
selected by our criteria. The continuous blue line is the 
sum of the probability distribution function in redshift for 
each BM/BX/LBG galaxy, derived by our photometric 
redshift code. It is in agreement with the distribution us- 
ing the best estimate values for the photometric redshifts 
(solid histogram). 



starbursts at intermediate/high redshifts. The BzK crite- 
rion is highly efficient in the redshift range 1.4 < z < 2.5, 
but when galaxies start to become faint, it is critical to 
distinguish between star-forming and evolved galaxies, re- 
sulting in an underestimation of the passively evolving 
population. 

3.4. Intrinsic properties 

To overcome the limitations of LBG, DRG and BzK selec- 
tion criteria, which are based simply on observed colours, 
we use another approach based on spectral fitting tech- 
nique. 

The spectral fitting technique adopted here is 
the same that has been developed in previous 
papers (jFontana et al. 2003). IFontana et al. 2004], 
IGrazian et al. 2006al and [Fontana et al. 2006), and 
similar to those adopted by other groups in the literature 
(e.g. IDickinson et al. 20031 |Drory et al. 2004[ ). Briefly, 
it is based on the comparison between the observed 
multicolour distribution of each object and a set of 
templates, computed with standard spectral synthesis 
models (Bruzual & Chariot 2003 in our case), and chosen 
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to broadly encompass the variety of star-formation 
histories, metallicities and extinction of real galaxies. To 
compare with previous works, we have used the Salpeter 
IMF, ranging over a set of metallicities (from Z = 0.02.Zq 
to Z = 2.5Z G ) and dust extinction (0 < E(B — V) < 1.1, 
with a Calzetti or a Small Magellanic Cloud extinction 
curve). Details are given in Table 1 of lFontana et al. "2 004 
For each model of this grid, we have computed the ex- 
pected magnitudes in our filter set, and found the 
best fitting template with a standard \ 2 normalisation. 
The stellar mass and other best-fit parameters of the 
galaxy, like SFR, age, r (the star formation c-folding 
timescale), metallicity and dust extinction, are fitted 
simultaneously to the actual SED of the observed galaxy. 
The redshift of each galaxy is fixed during the fitting 
process to the spectroscopic or photometric redshift 
derived in IGrazian et al. 2006al Clearly, this approach 
requires the availability of spectroscopic or photometric 
redshifts of good quality, in order to derive precisely rest 
frame properties for these galaxies. 

We define passively evolving galaxies according to the 
physical criterion age/r > 4. This quantity is in prac- 
tice the inverse of the Scalo parameter and a ratio of 4 is 
chosen as an arbitrary value to distinguish galaxies with 
pure evolved stellar populations from galaxies with recent 
episodes of star- formation. Moreover, an age/r — 4 corre- 
sponds to a residual SFR 2% of the initial SFR, for an ex- 
ponential star formation history, as adopted in this paper. 
We have checked that around this typical value galaxies at 
2 < z < 3 show a prominent 4000 A and Balmer break, the 
typical signatures of an emerging old stellar population. 

The two parameters used in this context, the age/r ra- 
tio and the stellar mass, are subject to uncertainties and 
biases related to the synthetic libraries used to carry out 
the fitting of the galaxy SEDs, as we discuss here. In ad- 
dition, the stellar mass generally turns out to be the least 
sensitive to variations in input model assumptions, and the 
extension of the SEDs to mid-IR wavelengths with IRAC 
tends to reduce the formal uncertainties on the derived 
stellar masses, as shown in IFontana et al. 20061 

The uncertainties in the ratio age/r are derived as 
follows: we compute the la confidence level on each esti- 
mate of age and r parameters by scanning the \ 2 levels, 
allowing the redshift to change in case of objects with pho- 
tometric redshifts. The typical uncertainty, shown in Fig. 
01 for the age/r parameter is of the order of 50%, while 
the uncertainty for the mass estimate within the GOODS- 
MUSIC sample is 40%. Though the age and r parameters 
are individually poorly constrained and highly degenerate, 
the ratio age/r is comparatively better constrained. 

Using the criterion based on the age/r ratio on our Ks 
selected sample, we identify 130 passively evolving galax- 
ies [age/r > 4) in the redshift range 1.4 < z < 2.5. 
Analogously, we find 508 young star forming galaxies (de- 
fined by the criterion age/r < 4) in the same redshift 
range. 

Fig. O compares these galaxies selected through the 
age/r index with those identified by the BzK-PE and 
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Fig. 4. The relative uncertainties of the age/r ratio. The 
dashed histogram refers to galaxies with known spectro- 
scopic redshift, whose age/r parameter has a relative error 
of 40%. The typical uncertainty of this ratio for galax- 
ies with a photometric redshift determination is instead 
~ 50%. 



BzK-SF criteria. It is clear from this plot that the BzK 
and age/r criteria are largely consistent. In particular the 
BzK-SF criterion is very effective in recovering the star- 
burst phase of galaxies at z ~ 2 since it is able to se- 
lect the large majority (94%) of galaxies with age/r < 4. 
On the other hand the BzK-PE criterion is more incom- 
plete (34%), since it misses a great fraction of passively 
evolving galaxies, having age/r > 4, especially at faint 
Ks magnitude limits and for objects with extreme z — Ks 
colours, as shown in the previous section. In this case, the 
DRG criterion turns out to have enhanced completeness 
in recovering passively evolving galaxies, but also strong 
contamination from dusty starbursts at 1 < z < 2. As 
a further validation of our approach based on the age/r 
ratio, we have selected all passively evolving galaxies with 
a <?e/r > 4 and 1.4 < z p hot < 2.5 in the BzK-SF locus and 
found one galaxy (ID=9853 in GOODS-MUSIC, ID=139 
in K20) with a spectroscopic redshift of z = 1.553 from 



the K20 survey flMignoli et al. 2005[ ) without [Oil] line in 
emission and with Mgll line in absorption; it is thus con- 
sistent with being an evolved galaxy. Other galaxies of 
this type are fainter, but the spectroscopic identification 
is at the reach of deep or ultradeep spectroscopy carried 
out by the GO ODS (IVanzella et al. 2006|) and GMASS 
(jCimatti 20061) teams in the CDFS field. 

We investigate also the distribution of age/r for differ- 
ent galaxy types. The distribution of the age/r ratio cal- 
culated using the probability distribution function of this 
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Fig. 5. The B—z vs z—Ks colour for ifs-selected galaxies 
in the GOODS-MUSIC sample, as in Fig. H Green filled 
squares show the distribution of all ifs-selected galaxies 
with redshift (spectroscopic or photometric) between 1.4 
and 2.5, with young stellar populations (age/r < 4). Cyan 
circles are galaxies with 1.4 < z < 2.5 and age/r > 4, a 
more efficient way to isolate passively evolving galaxies 
than the original BzK-PE criterion. 



parameter for each galaxy, is consistent within the uncer- 
tainties with the histogram built using the best fit values 
for each object. This shows that our approach, based on 
the best fit quantities for age/r are equivalent to results 
that one would obtain with a Monte Carlo approach. Fig. [6] 
confirms that BzK-SF and DRG galaxies are mixed popu- 
lations, and the distribution of the age/r parameters is an 
extended function, in which the relatively young objects 
dominate, especially for the BzK-SF class. For DRGs the 
distribution is less steep, with 110 objects with age/r < 4 
and 69 objects with age/r > 4. The BM/BX/LBGs are 
almost all at age/r < 4, confirming that they are pre- 
dominantly young star forming galaxies, while the BzK- 
PE class, despite the low statistic, shows a bimodal dis- 
tribution, with equal populations of young starbursts and 
passively evolving objects. These facts strengthen our pre- 
vious conclusions. 

Note that from these figures ([5] and [5]) we infer a poste- 
riori that the age/r — 4 value is a good separator since an 
apparent bimodality is evident in the age/r distribution. 
In particular, as shown in Fig. \6\ the age/r distribution 
for the Ks-selected galaxies flattens at age/r > 3.5 — 4.0, 
indicating that a population of passively evolving galaxies 
can be isolated with this proposed criterion. 



Fig. 6. Upper panel: The distribution of the age/r param- 
eter for the total Ks-selected sample, BM/BX/LBG-Z and 
DRG. The continuous solid curve shows the distribution 
of the age/r ratio calculated using for each galaxy the 
probability distribution function of this parameter. The 
curve is consistent within the uncertainties with the his- 
togram built using the best fit values for each object. This 
shows that our approach, based on the best fit quantities 
for age/r is equivalent to the result that one would obtain 
with a Monte Carlo approach. Lower panel: The same of 
the upper panel for the total Ks-selected sample, BzK-SF 
and BzK-PE galaxies. 



4. The stellar mass and mass density of high 
redshift galaxies 

The contribution of massive galaxies to the stellar 
mass function at high redshifts has been already dis- 
cussed by Fontana et al. 20061 |Caputi et al. 2006 



and 



IMcLure et al. 20061 A study of IFranceschini et al. 20061 
based on detailed morphological analysis with HST of a 
sample of galaxies detected at 3.6 fim in the GOODS- 
South, indicates that the high mass tail of the stellar mass 
function at z ~ 1 — 2 is due mainly to evolved galaxies. 
It is thus interesting to derive the relative contribution of 
DRGs and BzK-PEs to the stellar mass density (SMD) at 
high redshift. To provide a clear comparison, we compute 
the SMD also for LBGs and BzK-SF galaxies, deriving the 
contribution of star forming objects to the universal SMD 
at high-z. Finally, we also derive the mass distribution and 
mass density for age/r < 4 and age/r > 4 galaxy sam- 
ples, and provide a global picture on the SMD for various 
galaxy types at z > 1. 

The method that we applied to estimate the galaxy 
stellar masses is described in Section 3.4. A cautionary 
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note on the stellar mass derivation is necessary here, 
to remind that it always depends on the assumed star- 
formation history. As discussed by many authors (e.g. 
IFontana et aI_2004|and |Shapley et al. 2005[ ), the assump- 
tion of an exponential decreasing SFR for all galaxies may 
provide biases in the estimate of the stellar mass. In par- 
ticular, it may lead to an underestimate of the contribu- 
tion of older stellar population in the case of an actively 
star-forming galaxy, whose luminosity is dominated by 
the glare of young stars. The advantage of the exponential 
approach is that it can be compared with many previous 
estimates. In Fonta na et al. 20041 we have shown that it 
may lead to a maximal underestimate of the stellar mass 
by a factor of 2, in the most extreme cases. Another ap- 
proach is the so-called maximally old model, in which the 
near- and medium-IR SED is fitted with a burst of star 
formation at age ~ 0, while the UV part of the spectrum 
is reproduced by a more recent burst. This method recov- 
ers a larger stellar mass than that evaluated by the best fit 
approach, but is based on extreme assumptions, as shown 
in IFontana et al. 2004[ and it is only useful to provide an 
upper limit to the stellar mass assembled in a galaxy. 
We finally note that recent papers (|Rettura et al. 2 006: 
lErb et al. 2006P show that the dynamical mass derived 
from spectroscopy is consistent with the stellar mass de- 
rived by the best fit approach, reinforcing the robustness 
of our method. 



4.1. The observed mass distribution of high redshift 
galaxies 

We first present the observed distribution of stellar 
masses. Fig. [Jj (left) represents the distribution of stel- 
lar mass (note the logarithmic scale in the y-axis) for 
BM/BX/LBGs-Z, BzKs and DRGs in the redshift interval 
1.4 < z < 2.5, while Fig. [7] (right) refers to the redshift 
interval 2.5 < z < 3.6. Since the magnitude limit of the 
GOODS-MUSIC sample is not unique both in z- and in 
Ks-band, the resulting mass distributions are weighted for 
the different magnitude limits in the various sub-areas (see 
Section 2 of this paper and also IGrazian et a l. 2006a) . 

It must be emphasised that the observed distributions 
do not arise from a mass-limited sample, but rather from 
magnitude-limited samples (in the z or in the Ks-band), 
and are therefore to be critically used. At high masses, all 
the samples are complete (in stellar mass) : all the conclu- 
sions drawn on such samples at high masses are therefore 
robust. At low masses, two different effects are in place. 
First, the z-selected sample extends to fainter flux levels, 
and hence lower masses. In addition, even within the Ks- 
selected samples (DRGs and BzKs), star-forming spectral 
types extends to lower masses because their typical M/L 
ratio is lower than those of passively evolving objectfl 



We remind that the stellar mass of each galaxy is derived 
through the mass-to-light ratio fitted by the spectral synthesis 
model. Since evolved galaxies and star forming galaxies are 
characterised by distinct M/L ratios, a given cut in luminosity 



As described at length in IFontana et al. 20061 a safe mass 
threshold, i.e. the mass over which the samples are def- 
initely (~ 100%) complete in stellar mass, may be com- 
puted at each redshift using a maximally old, single burst 
model. When applied to the present if s-selected sample, 
it turns out to be log(M*/M ) ~ 10.6 at z ~ 2 and 
log(M*/M ) ~ 10.85 at z ~ 3. For the BM/BX/LBGs, 
selected from the z band, we do not apply any cut in stel- 
lar mass, since the relation between rest frame UV light 
(redshifted in the optical bands) and the stellar mass is 
too scattered at these redshifts for this type of galaxies. 

The observed mass distributions are also sensitive to 
the redshift selection function of different galaxy types. 
This selection function is not known a priori, but from 
the observed redshift distributions, which are not uniform 
over 1 < z < 4, as shown in Fig. [31 it is clear that the 
selection functions of BzKs, DRGs and BM/BX/LBGs are 
quite different. 

We therefore emphasise that the results of this sub- 
section are intended to provide a comparison of the con- 
tribution of different samples of galaxies as provided by 
current methods and surveys. A comparison of their over- 
all impact on the global stellar mass density is discussed 
instead in the following section. 

With all these caveats in mind, we can compare the 
observed distributions of stellar masses. At 1.4 < z < 2.5 
(Fig. O left), the fraction of galaxies at the high-mass tail 
(M > lO n M ) is dominated by DRGs and by BzK galax- 
ies (both SF and PE), while LBGs dominate the distribu- 
tion at lower masses. Current surveys detect an extended 
tail of low-mass star-forming galaxies (LBGs and BzK- 
SF) that is missed if the mass budget at high-z is measured 
only using DRGs. Such galaxies are much more numerous 
than the higher mass, redder galaxies, and their contri- 
bution to the mass budget is not negligible, as we shall 
discuss in the following. At higher redshifts (Fig. [71 right), 
the mass distributions of both DRGs and BX/LBGs re- 
semble those at lower redshift, especially in the massive 
tail. 

We also investigate the stellar mass distribution ac- 
cording to the age )t parameter: in Fig. [8] we plot the 
distribution of stellar mass for age/r > 4 and age/r < 4 
galaxies in the redshift range 1.4 < z < 2.5 (left panel) 
and 2.5 < z < 3.6 (right panel). 

Fig. [8] shows that at M > 10 n M Q and at z ~ 2 
the passively evolving galaxies are the dominant pop- 
ulation, while at z ~ 3 the high mass tail is mainly 
due to star forming galaxies. We are thus witnessing the 
epoch of mass assembly ("upsizing") of passively evolv- 
ing galaxies, confirming the strong density evolution of 
the Stellar Mass Function for such galaxies at z > 2, as 
shown in Drory et al. 2005 [JYanccschini et al. 2006 and 
IFontana et al. 20061 On the other hand, the high mass tail 
of star forming galaxies (age/r < 4) shows no evolution 
from z ~ 2 to z ~ 3, while the low mass distribution ex- 



translates in different cuts in mass for evolved and young stellar 
populations 
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Fig. 7. The distribution of stellar mass (in M Q and logarithmic scale) for DRGs, BzKs (SF and PE) and BM/BX/LBGs 
in the redshift range 1.4 < z < 2.5 (left panel) and 2.5 < z < 3.6 (right panel). Error bars represent uncertainties at 
68% confidence limit. 



periences a strong negative evolution at higher rcdshifts, 
yet an other indication of the so-called "downsizing" . 

We finally note that the mass distribution of passively 
evolving galaxies extends somewhat below the mass cut 
of the BzK-PE. The fall down of the mass distribution 
for BzK-PE is indeed due to the selection effect already 
described in Sect. 3.5, by which the BzK criterion turns 
out to be incomplete in the selection of passively evolving 
galaxies at faint Ks magnitudes. 

4.2. The stellar mass density of high redshift galaxies 

The Stellar Mass Density analysis for high-z galaxies is a 
more appropriate description than the mass distribution 
presented in the previous paragraph, since it takes into 
full account volume effects. In addition, and even more 
important, it allows to compare the contribution of the 
different classes of galaxies, as selected with the various 
methods described here, with the total stellar mass den- 
sity, as obtained by integrating the stellar mass function 
at the corresponding redshifts. 

The contribution to the total observed SMD at 
2 < z < 3 and M > 10 n M Q has been de- 
rived by Ivan Dokkum et al. 20061 and is 77% and 17% 
for DRGs and LBGs, respectively. In a similar way, 
IRudnick et al. 20061 measured the SMD of DRGs at z ~ 
2.8 over different lines of sight and found that these galax- 
ies contribute to 64% of the observed SMD at this redshift. 

We compute the contribution to the SMD for LBGs 
and DRGs in the same mass and redshift range of 
Ivan Dokkum et al. 20061 and we find that BM/BX/LBGs- 



Z recover 9% of the total observed SMD, while DRGs 
make up 64% of the observed SMD. This result agrees 
quite well with that of Ivan Dokkum et al. 20061 and 
IRudnick et al. 20061 

We can expand this analysis using our deeper GOODS- 
MUSIC catalog: we derive the mass density p* M for DRGs, 
BzKs and BM/BX/LBGs separately as a function of red- 
shift, and compare it with the observed p* M for the Ks 
selected sample and with the total value obtained in- 
tegrating the fitted mass function over the mass range 
10 8 - 10 12 M Q (jFontana et al. 2006|) . We then repeat the 
same analysis on the age/r selected samples. 

We have divided the DRGs and BM/BX/LBGs sam- 
ples in three redshift intervals, following the BM, BX and 
LBG redshift distributions. The BzK contribution to the 
stellar mass density is computed only in the redshift in- 
terval 1.4 < z < 2.5 where the BzK criterion is effective. 
The age/r selected samples are arbitrarily divided into 
two redshift intervals, 1.4 < z < 2.5 and 2.5 < z < 3.6; 
the first redshift bin is chosen so that direct comparison 
with the BzK criterion is possible, while the second bin is 
almost identical to the high redshift DRG and LBG red- 
shift distributions as in Fig.[7]and[8l Fig.[9]shows the SMD 
for DRGs, LBGs and BzK galaxies at different redshifts, 
while Tab. [5] contains the SMDs for these galaxy types. 
The error bars of the SMD have been computed with a 
full Monte Carlo simulation where we take into account 
the redshift probability distribution for each galaxy in the 
sample. 

These numbers are to be compared with the to- 
tal mass density obtained by integrating between 10 s 
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Fig. 8. The distribution of stellar mass (in M and logarithmic scale) for age/r > 4 and age/r < 4 galaxies in the 
redshift range 1.4 < z < 2.5 (left panel) and 2.5 < z < 3.6 (right panel). Error bars represent uncertainties at 68% 
confidence limit. 



and 1O 12 M the global Galaxy Stellar Mass Function 
(GSMF) observed in the same GOODS-MUSIC sample 
(Fontan a et al. 2006|) . Clearly, these fractions are rela- 
tively small since the lower integration limit for the total 
SMD is far below the observational limits for both the Ks 
and the z-selected samples. Finally we repeat the same 
calculations for M > 10 n M©. 

The value of the total (integrated) stellar mass den- 
sity depends on the slope of the GSMF at the low mass 
side, which is presently not known at these redshifts. The 
Schechtcr slope index a is observed to remain considerably 
flat, changing from a = —1.18 at z = to a ~ —1.3 at 
z ~ 1.3 (IFontana et al. 2006)) . As a first estimate, we shall 
assume that this trend of slowly decreasing a will con- 
tinue up to z ~ 3 (where a would be —1.54), as implicit 
in the Schechter parametric representation of the GSMF 
provided in Fontana et al. 20061 We shall discuss in the 
following the impact of releasing this assumption on the 
slope of the GSMF. 

In the first redshift bin, 1.4 < z < 2.2, DRGs recover 
21% of the total mass density, similar to the contribution 
of BM-Z galaxies, that is 14% (see Tab.[2]). BzK-SF proves 
to be a more efficient criterion recovering 75% of the to- 
tal SMD, while BzK-PE provide around 19% of the total 
mass density. Consistently with our previous discussion of 
the relative efficiency of the selection criteria, the BzK cri- 
teria turn out to be more efficient than BM-Z or DRGs for 
deriving the total SMD. Moreover, BzK-SF galaxies con- 
tribute to 75% of the total SMD, indicating that at these 
redshift a considerable fraction of the mass budget can be 
found in star forming galaxies, possibly obscured by dust. 



In the higher redshift bins, BX/LBG-Z galaxies con- 
tribute to 58% and 48% of the total SMD at 2.2 < z < 2.8 
and 2.8 < z < 3.7, respectively, while DRGs recover 
37 and 28% in the same redshift intervals. This im- 
plies that BX/LBG-Z galaxies cannot be neglected when 
the total SMD is investigated. Of course, limiting the 
analysis to the high mass tail of the galaxy mass func- 
tion, i.e. masses larger than 10 11 M Q , the contribution of 
DRGs is predominant, as we can see in the lower part of 
Tabj2]and in agreement with Ivan Dokkum et al. 20061 and 
IRudnick et al. 20061 

To further check these results we have derived the con- 
tribution of passively evolving and star forming galaxies 
selected according to the age/r parameter: star forming 
galaxies provide a considerable fraction of the total mass 
budget, in particular at redshift ~3. The contribution of 
passively evolving galaxies (age/r > 4) at 1.4 < z < 2.5 is 
32%. At higher redshift the contribution of such galaxies 
becomes even lower, only 13%. 

These results are in practice obtained by integrating 
the observed mass distributions of Fig. [7] and Fig. [8] (with 
the proper volume element) and comparing them with the 
integrated GSMF. As mentioned before, the observed dis- 
tributions of Fig. [7] and Fig. [5] are biased by selection ef- 
fects against low mass, red galaxies, while the total SMD 
depends on the assumed slope of the GSMF. Because of 
these systematics, it is possible that the observed distri- 
butions of Fig. [7] and Fig. [5] are only the results of a selec- 
tion effect, and that red (or passively evolving galaxies) 
keep dominating the stellar mass budget even at lower 
masses, overpredicting the total SMD. This would imply 
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that the slope of the GSMF becomes significantly steeper 
than what assumed here at high redshift. In the following 
subsection, we try using the same GOODS-MUSIC data 
as a first attempt in this direction. 

4.3. Looking for faint DRGs in the Ks-band 

As discussed above, the DRGs and in general all red galax- 
ies are limited to high masses simply due to their typi- 
cally higher M/L: pushing our analysis to one magnitude 
deeper than the actual magnitude limits, it is possible to 
investigate the existence of a population of faint/low-mass 
DRGs, whose contribution to the SMD could overcome 
that of LBGs. 

To push the present sample toward lower mass galax- 
ies, we derive the stellar mass density shifting by one 
magnitude both the actual magnitude limits of the sur- 
vey, in the z as well as in the Ks bands, (corresponding 
to 90% completeness all over the GOODS-South field, as 
described in Section 2). In this case, the distribution is 
only a lower limit to the number density of galaxies at 
low masses, while it is well constrained at the high mass 
tail. If we compute the SMD for LBGs and DRGs using a 
mass limit ~2.5 times lower (corresponding to one magni- 
tude deeper in the z and Ks bands, or alternatively to 75% 
completeness) than the conservative one, we derive a SMD 
higher both for LBGs and DRGs, as we expect. The differ- 
ential increase however is not constant, as shown in Fig. 
(upper panel), but DRGs gain more (~ 1.6 times) than 
LBGs in terms of SMD, especially at z ~ 3. This provides 
an indication for the existence of a numerous population 
of faint and less massive DRGs than those observed today, 
although it is still unclear whether they could overcome 
in number and mass density the contribution of LBGs at 
high redshift. 

The mass distribution of DRGs thus seems to be con- 
tinuous and not limited to the high mass tail. A survey 
much deeper than GOODS, like the ultra-deep FIRES 
Ks-band observation of the HUDF (jBeckwith et al. 2003| 
Labbe et al., in prep.), will probably be able to constrain 
the number density of low-mass DRGs in this redshift 
range. 

5. Discussion and Conclusions 

In this work we have used the GOODS-MUSIC sample to 
compare the colour criteria commonly used for selecting 
high redshift galaxies, and to estimate their contribution 
to the universal stellar mass budget. 

The GOODS-MUSIC sample is particularly suited 
for this exercise, since it provides a reasonable statistics 
(nearly 3000 galaxies in the ifs-selected sample and more 
than 9000 in the z-selected one), a 14 bands wavelength 
coverage that allows a direct application of most colour se- 
lection criteria and well tested, accurate photometric red- 
shifts. Most important, it includes the 3.6 — 8fim Spitzer 
observations of the complete data set, that are necessary 
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Fig. 9. Lower panel: the stellar mass density of the ob- 
served Ks-selected sample (circles), DRGs (filled trian- 
gles) and LBGs (filled squares) for the redshift bins where 
the BM, BX and LBG selections are efficient. DRGs sam- 
ple a considerable fraction of the stellar mass density of 
the Universe at z ~ 2.5, but LBGs are complementary to 
DRGs in recovering the total p* M . The BzK-SF criterion 
(void triangles) is very efficient in the lower redshift bin 
1.4 < z < 2.5, surpassing also the contribution of the DRG 
sample. The contribution of BzK-PE galaxies (asterisk) is 
slightly reduced, due to incompleteness of this selection 
criterion. Contrarily, the selection age/r > 4 (crosses) re- 
covers a higher stellar mass density, which is the contri- 
bution of true passively evolving galaxies at z ~ 2. The 
dominant sample at z > 1 is however that selected by 
a 9 e / T < 4 (stars). The continuous line shows the mass 
density obtained integrating the total galaxy stellar mass 
function between 10 8 and 10 12 M©. Upper panel: the mag- 
nitude limits both in the z and in the Ks bands are then 
lowered down of one magnitude in order to try to derive 
the contribution of faint galaxies to the mass density at 
high-z (pdeep)- At z ~ 3 the increase on SMD for faint 
DRGs is larger by a factor ~ 1.6 than the increase due to 
faint LBGs. 



to provide reliable estimates of the stellar mass for galax- 
ies at relatively high redshift (Font ana et al. 2006). 

We have initially selected samples of BM/BX/LBGs, 
DRGs and BzK galaxies, discussing the overlap between 
the various sets and the limitations of each criterion. As 
expected from previous works, the selection of galaxies ac- 
cording to the BM/BX/LBG criteria is sensitive to mod- 
erately obscured star-forming galaxies, missing the dusty 
starburst objects. DRGs, instead, are less sensitive to dust 
obscuration effects, but comprise a mix of two populations, 
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Table 2. Stellar mass density for DRGs, BM/BX/LBGs 
and BzK galaxies 



Sample 


a ~ 1.8 


z ~ 2.5 


z ~ 3.3 


Log p\[ 


TOT 


7.947 


7.681 


7.380 


OBSERVED-Ks 


7.924 


7.485 


7.118 


DRG 


7.275 


7.248 


6.823 


BM/BX/LBG-Z 


7.096 


7.440 


7.064 


BzK-SF 


7.774 






BzK-PE 


7.177 






age It > 4 


7.402 




6.601 


age )t < 4 


7.728 




7.096 


Type/TOT 


DRG 


0.21 


0.37 


0.28 


BM/BX/LBG-Z 


0.14 


0.58 


0.48 


BzK-SF 


0.75 






BzK-PE 


0.19 






age/T > 4 


0.32 




0.13 


age/T < 4 


0.68 




0.42 


Type/TOT(M > 10 lx M e ) 


DRG 


0.30 


0.77 


0.79 


BM/BX/LBG-Z 


0.01 


0.09 


0.10 


BzK-SF 


0.57 






BzK-PE 


0.37 






age/T > 4 


0.55 




0.25 


age/T < 4 


0.45 




0.75 



p* M is the stellar mass density in MqMpc 3 . The three red- 
shift intervals, 1.4 < z < 2.2, 2.2 < z < 2.8 and 2.8 < z < 3.7, 
correspond to the BM, BX and LBG selection criteria, respec- 
tively. p* M (TOT) is the mass density obtained integrating the 
stellar mass function between 10 s and 10 12 Mq. The BzK cri- 
terion is effective only in the redshift range 1.4 < z < 2.5, and 
for the higher redshift bins it is not possible to compute its 
stellar mass density. The SMD for galaxies divided according 
to their age/T parameter is computed in the redshift intervals 
1-4 < z < 2.5 and 2.5 < z < 3.6. The total SMD in these red- 
shift ranges are 10 7,898 and 10 7 ' 471 , and these two quantities 
are used to compute the relative contribution of BzK galaxies 
and of objects selected by their age/T ratio. 

the old/evolved galaxies and the dusty starbursts at inter- 
mediate/high redshifts. The BzK criterion is highly effi- 
cient in the redshift range 1.4 < z < 2.5, but when galaxies 
start to become faint in the Ks band and red in the z — Ks 
colour, it is difficult to distinguish between star- forming 
and evolved galaxies, resulting in an underestimation of 
the passively evolving population. 

To better separate actively star-forming galaxies from 
passively evolving ones, we have then applied a physical 
criterion, based on the ratio between the age of the stel- 
lar component and the star-forming timescale, both de- 
rived from the best-fitting models applied to the full SED 
of each galaxy. We show that the ratio age/r is a rela- 
tively well constrained parameter, and we adopt a thresh- 
old age/r > 4 to separate passively evolving galaxies from 
star forming ones. 

We have then analysed the observed distribution of 
stellar masses associated to each selection criterion. At 
1.4 < z < 2.5, the fraction of galaxies in the high-mass 



tail (M > 10 n Af Q ) is dominated by DRGs and by BzK 
galaxies (both SF and PE), while BM/BX and BzK-SF 
dominate the distribution at lower masses. The same kind 
of bimodality is apparently in place at 2.5 < z < 3.6, 
where DRGs make the high mass tail of the distribution 
and BX/LBGs the low mass side. Clearly, the lack of red 
galaxies with low stellar mass results from the selection 
criteria adopted: since the z-selected sample extends to 
fainter fluxes than the ifs-selected one, and given that 
the M/L ratio of red galaxies is larger than that of blue 
galaxies, current samples are biased against the detection 
of low mass red galaxies. 

If we look at the intrinsic properties, it is remarkable 
that at z ~ 2 passively evolving galaxies, selected with 
age/r > 4, exist in a large number, such that they are 
slightly dominant at the highest masses (M > IO^Mq) 
over the star-forming population. However, they show a 
significant negative evolution with redshift, with a more 
modest contribution to the high mass tail at z ~ 3. This 
transition suggests that the epoch from z ~ 2 to 3 is 
crucial in the assembly of passively evolving galaxies. 

We have then derived the contribution of each se- 
lected sample to the Stellar Mass Density (SMD) at var- 
ious redshifts. This is compared with the total SMD, 
as obtained by integrating the Galaxy Stellar Mass 
Function of lFontana et al. 20061 from 10 8 Af o to 1O 12 M . 
Such a comparison is based on the assumption that 
the slope of the Galaxy Stellar Mass Function at high 
redshift remains relatively flat, as inferred from an 
the extrapolation of its evolution at lower redshifts 
(jFontana et al. 2006|) . Under this assumption, the over- 
all contribution of BM/BX/LBGs to the total SMD from 
z ~ 1.5 to z ~ 3.5 is comparable to or even higher than 
that of DRGs, since they outnumber the DRGs at low 
masses. Clearly, if only the high mass tail is considered, 
DRGs are the dominant population. The importance of 
star forming galaxies is even more noticeable if one con- 
siders the BzK-SF sample, which is dominating the SMD 
at 1.4 < z < 2.5, since it recovers both moderately ob- 
scured and dusty young galaxies. 

If this picture is correct, it presents a similarity and 
a difference with the present-day Universe. The similar- 
ity is that a bimodality seems to exist, with star-forming 
galaxies dominating the low mass population and pas- 
sively evolving galaxies dominating the high mass tail. At 
variance with the local Universe, however, the integrated 
contribution of passively-evolving galaxies is lower than 
that of star-forming ones. These two clues provide a dif- 
ferent view of the so called "downsizing" scenario. 

However, there is the possibility that a large popu- 
lation of low-mass, intrinsically red galaxies exists, still 
undetected in current K-selected surveys. Such a popu- 
lation could in principle produce a Galaxy Stellar Mass 
Function much steeper than assumed here, and would 
significantly change the relative contribution of different 
galaxy types to the SMD. Detecting the existence of such a 
population, and constraining its nature (passively evolving 
vs. star-forming) would be very important to understand 
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the physical processes in high redshift galaxies, since in 
theoretical scenarios feedback effects and star-formation 
histories are strongly dependent on the halo mass. Only 
much deeper IR-selected surveys will be able to discover 
whether the contribution of DRG and other red galaxies 
to the SMD at z > 2 is actually higher than what ob- 
served with present data: DRGs indeed could eventually 
overcome the BM/BX/LBG SMD when their contribu- 
tion is integrated till the very low mass tail of the mass 
distribution. 
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